Effective Mean Field Approach to Kinetic Monte Carlo Simulations in Limit Cycle 
Dynamics with Reactive and Diffusive Rewiring 
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The dynamics of complex reactive schemes is known to deviate from the Mean Field (MF) theory 
when restricted on low dimensional spatial supports. This failure has been attributed to the limited 
number of species- neighbours which are available for interactions. In the current study, we introduce 
effective reactive parameters, which depend on the type of the spatial support and which allow for an 
effective MF description. As working example the Lattice Limit Cycle dynamics is used, restricted 
on a 2D square lattice with nearest neighbour interactions. We show that the MF steady state 
results are recovered when the kinetic rates are replaced with their effective values. The same 
conclusion holds when reactive stochastic rewiring is introduced in the system via long distance 
reactive coupling. Instead, when the stochastic coupling becomes diffusive the effective parameters 
no longer predict the steady state. This is attributed to the diffusion process which is an additional 
factor introduced into the dynamics and is not accounted for, in the kinetic MF scheme. 

PACS numbers: 



I. INTRODUCTION 

The theory of reactive dynamics has for a long time 
mainly been restricted to studies of the macroscopic, phe- 
nomenological, Mean Field (MF) equations. As a result, 
effects such as local interactions, spatial restrictions, de- 
fects, local stochastic effects, etc were often ignored or 
added ad hoc. In the recent years, with the development 
of the Kinetic (or Dynamic) Monte Carlo Methods, it is 
possible to include in detail these factors and to follow 
the system as it evolves dynamically from one state to 
another [ll-0|- By generating a dynamical state-to-statc 
trajectory it is possible to explore the entire state space 
as the system is directed towards the steady state, while 
the dynamics and the steady state crucially depend on 
the dynamical history. 

In models of chemical catalytic dynamics S SJl, @- 
[Toj , ecological models flll - [l3| , epidemiology [MtIII, the 
existence of a spatial support is crucial and may modify 
considerably the MF approximation. In addition, diffu- 
sion of species may often modify the processes [ill, [111 ■ 
In all these systems the supports present certain degrees 
of complexity. For example, in ecological systems the dif- 
ferent species live and interact in natural environments 
that present differences in their structure from a sub-area 
to another [ll|, [13. As it is noted in [l3], in the chemi- 
cal, biological and ecological world, pattern formation is 
a very common phenomenon, that must be always taken 
into account in respective models, since they affect the 
dynamics between the species. 

In a classic work [^ the NO + CO reaction on Pt sur- 
faces was studied, without considering the corresponding 
spatial effects. The stability of the system was investi- 
gated and local kinetic oscillations were considered. In a 



recent investigation of the same system [l3|, it is noted 
that the dynamics leading to self-organization affects the 
" spatio-temporal organization" of the chemical catalytic 
reactions. 

In [ll| . synchronization among populations of differ- 
ent species which are related to each other as predators- 
consumers-vegetation, was studied. The dynamics of 
such systems is of the limit cycle type. It was found 
that there is a strong relation between the populations 
synchronization and the spatio-temporal characteristics 
of the system. This phenomenon takes also into account 
the possible "diffusive migration" of some species to a 
neighbouring or distant areas. 

In epidemiological models, the spatial characteristics 
of the systems are very important. In many studies, see 
e-g- [13 [i3|, it is shown that the structure of the social 
network affects the dynamics of a spreading disease [l3l , 
or that targeted immunization to the individuals-sites of 
the network that are best connected causes the localiza- 
tion of the disease [l^] . 

In order to study these spatiotemporal phenomena, a 
"center" dynamical system, the Lattice Lotka Voltera 
(LLV) system, was introduced on a 2-D square lattice 
[l8j . It was found that the conservative "center" dynam- 
ics reduces to local oscillations when the system is re- 
stricted on low dimensional supports. The attributes of 
the oscillations depend on the lattice size, the number of 
interacting neighbors and the general spatial restrictions. 
These oscillations are not of the limit cycle type, due to 
the nature of the interactions. Later on, long-distance 
diffusion [l^] was added in the LLV system and it was 
shown that in this case the out-of-phase local oscillations 
get phase synchronized and give global oscillations which 
are stable. After a critical point, a Hopf-like bifurcation 
happens and the system enters into limit cycle-like dy- 
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namics. 

A different model, involving a limit cycle MF dynamics 
on lattice was introduced in 2002, the Lattice Limit Cycle 
(LLC) model [lO, [Hj. For this model it was also shown 
that the behavior in the MF level and in the KMC simu- 
lations level are different. In the MF level the concentra- 
tions oscillate with constant amplitude, independent of 
the initial conditions (limit cycle dynamics). When the 
simulations are performed on a 2-D square lattice, local, 
out of phase oscillations, are observed between different- 
distant subregions of the lattice. Again, the spatial re- 
strictions of the system, the lattice size, etc, are strong 
parameters that affect the dynamical behavior of the con- 
centrations. 

In the current study, the addition of Long-Distancc Re- 
action (LDR) and Long-Distance Diffusion (LDD) pro- 
cesses is attempted on the LLC model. The aim of the 
study is to extend the previous works, by taking into 
account the spatial and stochastic coupling among the 
species involved in the reaction scheme and examine how 
the diffusion would alter the behaviour of these species 
and the dynamics of their concentrations. In addition, 
a posteriori effective parameters are introduced in the 
simulations which arc shown to restore the MF regime 
only in the case of LDR. In the case of LDD the MF re- 
sults different from KMC results even when the effective 
parameters are taken into account. 

In the next section the LLC model is presented where 
three different species interact in a predator-prey chain, 
and the corresponding MF equations are recapitulated. 
To implement the model. Kinetic Monte Carlo (KMC) 
simulations are performed on a square lattice and the re- 
sults are shown in Sec. |lll Differences between the MF 
and the KMC approaches are pointed out. In Sec. IIIII 
effective reactive parameters p'j^, P'3 ^-rc calculated a 
posteriori and the MF dynamics with effective parame- 
ters are compared successfully with the KMC results. In 
the next two sections two different mechanisms and their 
consequences are studied. In Sec. |lV]long distance reac- 
tions are added to the system, while in Sec. |V] long dis- 
tance diffusion is implemented. In both cases the species 
concentration is studied using the a posteriori effective 
values in the MF equations. The control parameters in 
these steps are the probability of long distance reaction 
and the probability of long distance diffusion. The con- 
cluding results as well as suggestions for future studies 
are presented in the Concluding section. 



II. THE LATTICE LIMIT CYCLE MODEL: 
MEAN FIELD AND KINETIC MONTE CARLO 
APPROACHES 

The LLC model is a model which describes the in- 
teractions among three completing species, in a chain 
of predator-prey interactions [13]. Each species acts as 
predator for one of the other two species. With Xi, X2 
the two interacting species (or particles) are denoted. 



whereas S is considered to be a fictitious -virtual species, 
representing the empty sites. In the model, the species 
live on a 2-D square lattice where each site can be occu- 
pied by a species Xi, or X2 or it can be empty. Single 
occupancy of all sites is only allowed at all times. The 
species S represents the empty sites - vacuum states of 
the lattice, which "interact" with occupied lattice sites 
as it will be seen in the next paragraph. 

Initially, the species are randomly distributed onto the 
lattice under given initial conditions. Each site can in- 
teract with its closest neighbours with given rates. The 
way that the species interact with each other is described 
in the following scheme: 

2X1 + 2X2 ^3X2 + 3 (la) 

Xi + S^ 2X1 (lb) 

X2 + S^2S (Ic) 

The reaction step [la] describes the main interaction 
between the species Xi and X2- It is a 4-th order non- 
linear reaction. In order for this reaction to take place 
two species Xi and two species X2 are needed. When 
these particles are found in the proximity of one another 
(closest neighbours) then with rate pi, one of the two 
Xi becomes an X2 while the second Xi desorbs (or dies) 
leaving the site empty (S). The reaction step I Ibl describes 
how the species Xi is born. The vacuum state interacts 
directly with particles (states) Xi and X2. In particular, 
when a particle Xi is found close to an empty site S*, 
then another Xi is born with reaction rate p2 and gets 
hosted on the empty site S which changes state from 
S ^ Xi. The reaction step [Tc] describes how the species 
X2 dies. When a particle X2 is found close to an empty 
site S then with rate pa, X2 dies leaving its site empty, 
X2 — > S . These interactions rates, which are given by 
the parameters pi, p2^ Pa, will be later translated into 
reaction probabilities when the simulation algorithm will 
be introduced. 



A. The Classical Mean Field Theory 

In this subsection we briefly recapitulate the main at- 
tributes of the classical MF theory describing the LLC 
scheme [13, [13 • The nonlinear kinetic equations for the 
species concentrations have the form: 

= -2^12:^X2 +^22^15 (2a) 
— = pix^x^ - p^X2S (2b) 

ds 2 2 /r> \ 

— =Pix^X2-p2Xis+p:iX2S. (2c) 
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where the small letters Xi, X2 and s represent the global 
MF particle concentrations. The space conservation con- 
dition 

Xi + X2 + s ~ const (3) 

is automatically satisfied. As usual, the constant is cho- 
sen equal to unity, leading to the interpretation of Xi, x-i 
and s as partial concentrations of particles Xi, X2 and 
empty sites. Using the conservation condition Eq. [31 it is 
straightforward to reduce the system by eliminating the 
s(= 1 — .Ti — X2) variable: 

^ = -2'Pyx\x\ ^ ■P2X\(\ - Xi ~ X2) 

dX2 2 2 \ / A\ 

— = piX^X2 ~ P3X2(l ~ Xi ~ X2j (4) 

This reduced system admits four steady state solutions, 
three of which are trivial and correspond to full occupa- 
tion of the lattice by one of the three species. The three 
trivial solutions can be represented as state vectors in 
the reduced {xi and X2) dimensions, namely Qi = (OiO) 
(empty lattice), Q2 = (0, 1) (lattice poisoned by X2) and 
Q3 = (1,0) (lattice poisoned by Xi). In addition the 
system [2] admits a fourth nontrivial solution with coexis- 
tence of all species in the steady state: 

\&^+\&^) 

where the constant K is only a function of the three 
reaction rates, 

K = Vl + (2p3+P2)V(27piP2P3). ^ 

Standard linear stability analysis indicates that the 
first three trivial fixed points are saddles, while the sta- 
bility of Q4 depends on the parameter values. For certain 
parameter values Q4 undergoes a supercritical Hopf bi- 
furcation, becomes an unstable focus and in its vicinity a 
stable limit cycle appears. In this regime the concentra- 
tions xi and X2 oscillate periodically, while the amplitude 
of the oscillations depends also on the system parameters. 

Note that although the MF description is appropri- 
ate for describing the main long time tendencies of the 
system, it is not a suitable approach when fiuctuations, 
spatial restrictions and stochastic effects are considered. 

B. The Kinetic Monte Carlo Approach 

When a detailed description of the system's dynam- 
ics is needed, where fluctuations, spatial restrictions and 
stochastic effects are considered, the ultimate descrip- 
tion is the probabilistic Master Equation approach. It 
describes the temporal evolution of the system from one 



state to another in detail. In particular, if the system is 
found in a finite number of states i = 1, 2, ...N then the 
probability P{j,t) to find the system in state j, at time 
t, is described by the Master Equation as: 

N N 

p{j, t)^Yl * - 1) - E ^^'^(j' * - 1) (6) 

k=l 1=1 

where the matrix element Wjk represents the transition 
probability from state j to state k. The matrix W with 
elements Wjk, j,k — 1 ■■■ N is called the transition prob- 
ability matrix. It is straight forward to see that the so- 
lution of the Master Equation is easy when the number 
of states is small, but in the case of systems with a 
large number of states, as the LLC lattice model, it be- 
comes impossible. An alternative method where the sys- 
tem samples its state space stochastically and thus for 
relatively large times gives a good approximating solu- 
tion to the Master Equation is the KMC method, which 
is used in the current study. 

The KMC method is a discrete time method starting 
from an initial random configurations of Xi, X2 and S 
particles on lattice and the update is random and sequen- 
tial, following the LLC reactive scheme ([T]). As substrate 
a 2D square lattice of size L x L is used and each particle 
interacts only locally with its nearest neighbours in the 
original scheme. Later, in sections IIVI and |Vl stochas- 
tic long distance couplings will be allowed with other, 
distant sites. Each lattice site, whose coordinates are de- 
noted as contains only one particle (Ai or X2) or 
is empty (S). Originally, the three species are randomly 
distributed on the lattice with given initial concentra- 
tions. Each Elementary Time Step (ETS) of the KMC 
algorithm contains four stages: 

1. One site, (i, j), is randomly chosen through out the 
lattice. 

2. If the site {i, j) contains an X2 particle and amongst 
the four nearest neighbours one particle of type S 
is found then the site changes its state from 
X2 to S with probability p^. This is the realisation 
of reaction Px)) . 

3. If the site contains an S particle and amongst 
the nearest neighbours a particle of type Xi is 
found then the site changes its state from S 
to Xi with probability p2. This is the realisation 
of reaction (jlbp . 

4. If the site (i,j) contains Ai, and its immediate 
neighbourhood contains one Ai and two A2 parti- 
cles, then the original and the neighbour Ai change 
simultaneously to S and A2, with probability pi. 
This is the lattice realisation of reaction ([Ta]) . 

5. If none of the above three steps is realised, the lat- 
tice remains unchanged. 

6. One ETS is completed. The algorithm returns to 
stage (1) for a new ETS to start. 



The indicative Monte Carlo time Step (MCS) consists of 
a number of elementary steps equal to the lattice size 
namely 1 MCS ^ ETS. In one MCS step each 
particle has reacted once, on average. The above KMC 
algorithm is a slight variance of the original one proposed 
in ref. [lO] , in order to establish consistency between the 
short distance interactions introduced in [20| and the long 
distance ones which will be introduced in the next two 
sections. 

As working parameter set the values pi = 0.9585, 
P2 — 0.016 and = 0.026 will be used which are lo- 
cated well inside the Hopf bifurcation region and giv e 
rise to periodic oscillations with large amplitudes [20|. 
In addition the working parameter values were chosen 
to be Pi < 1, so that they can directly be interpreted 
as reaction probabilities in the KMC scheme. In cases 
where the reaction rates are greater than unity they can 
be transformed into probabilities in two ways: a) ei- 
ther by dividing each one of them by the sum of all 
Pi (i.e. pi Pi/Y.Pj^ = 1:2,3) or b) by div- 
ing each Pi by the maximum of the reactive rates (i.e. 
Pi ~> Pi/p''"'"'^, i — 1, 2, 3). In both cases, this rate rescal- 
ing is equivalent to a time rescaling of the system j20j . 
Note that the large gap in the working set value of pi in 
relation to p2 and ps balances the highly selective, non- 
linear structure of reaction step (|lap . 

In Fig. [T]the average concentration of species Xi over 
all the lattice is depicted as a function of time for two dif- 
ferent lattice sizes, L = 2^ (red dashed line) and L = 2^" 
(blue solid line). As was also noted in [201 the restric- 
tion of the reactive dynamics on a support, the effects 
of spatial particle distribution and the stochastic noise 
induce the following modifications to the MF behaviour: 
a) the oscillations loose their regularity and become in- 
termittent (stochastic effects), b) the amplitude of the 
intermittent oscillations shrinks as the system size L in- 
creases, and this is attributed to the stochastic effects 
which randomise the phases of the local limit cycle os- 
cillators and c) a small shift is observed in the center of 
the KMC cycles with respect to the MF solution. For 
comparison the solid black line denotes the MF position 
of the center of the limit cycle. 




time (MCS) 



FIG. 1: Time series of average Xi density for different 
lattice sizes shows intermittent oscillations. The 
parameter values are: pi = 0.9585, p2 = 0.016 and 
Ps = 0.026. The solid black line denotes the MF steady 
state solution, Eq. [S] 



ing the KMC simulations. In many cases, where a parti- 
cle is selected for reaction with favorable rate, the event 
might not take place due to inappropriate environment 
(nearest neighbouring particles) or to stochastic choice. 
This way the original rates are modified and new, "ef- 
fective" rates govern the dynamics of the system. These 
rates are computed, a posteriori, within the algorithm 
provided in Sec. IIIBI 

During the application of the algorithm we introduce 
three event counters A^i, N2 and N3, which increase by 
one unit when the reactions [Tal [lb] and [Tc] take place, 
respectively. The "effective" parameters, pj, i ~ 1,2,3, 
are calculated as 1231: 
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III. CALCULATING THE EFFECTIVE 
PARAMETER VALUES 

As was shown in Sec. [Ill the stochastic and spatial 
effects modify the reactive dynamics with respect to the 
MF behaviour. Apart from the stochasticity in the ampli- 
tude of the oscillations, the position of the center, one of 
the most essential characteristics of the MF, is displaced. 
Because the position of the center is solely defined by the 
reactive parameters pi, p2 and ps, it is important to ex- 
plore whether the spatial and stochastic effects influence 
these parameters. 

To this purpose, we calculate, a posteriori, the "effec- 
tive" reactive parameters, which realistically occur dur- 



As an example, in Fig. [21 the effective values p'^ is plot- 
ted for various values of pi when the other parameters p2 
and P3 are kept constant. From the figure it is clear that 
as the reaction probability increases the effective reaction 
probability, which averages out all the local effects, also 
increases proportionally. 

Using the effective parameters p'^, the effective MF 
steady state values 2 and s' are calculated using Eq. 
[5l where the kinetic parameters are replaced with the 
effective ones. The results are plotted in [3| for various 
values of the original parameter pi and keeping the other 
two constant (p2 = 0.016 and pz = 0.026). In particu- 
lar, in Fig. [3|the three curves represent: a) the averaged 
values {xij taken directly from the KMC simulations, av- 
erage taken over 10 runs (black line with circles), b) the 
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FIG. 2: Average effective reaction probabilities p'l as 
functions of the imposed reactive reactive coupling pi . 

The other parameter values are: P2 = 0.016 and 
P3 = 0.026 and the system size is L = 2^. Averages are 
taken over 10 runs. The error bars on the p'l values are 
of the order of 4 x 10^^. 



MF calculated effective partial density x'^ of species Xi, 
using the effective probabilities p[ in the MF Eq. [Sj (red 
line with crosses) and c) the original MF partial den- 
sity xi, using the original probabilities pi in the MF Eq. 
[5] (green line with diamonds). The simple MF clearly 
underestimates the average partial population densities, 
while the effective MF gives a very close approximation 
to the KMC simulations. These results indicate that, as 
far as the steady state properties of the LLC model are 
concerned, the restriction of the system on the substrate 
together with the stochastic character of the reaction lead 
to new effective parameter values, while the MF steady 
state is achieved for these shifted parameter values. 



IV. KINETIC LATTICE MONTE CARLO WITH 
STOCHASTIC REACTIVE REWIRING 

As discussed earlier, stochastic rewiring may take place 
in networks where species interact with far away "neigh- 
bours" . As an example we give opinion exchange through 
the Internet or through the telephone network. Opinions 
are discussed and exchanged without actual displacement 
of the individuals. In such interaction environments the 
individuals are exposed both to the local environment 
(family, work) where the interactions are governed by 
rates pi , and to the long distance environment where the 
interactions arc governed by rates Pi.iong which, in gen- 
eral, are different. The long distance reactions introduce 
a type of reactive mixing in the system and drive the 
system towards its MF behaviour. In some cases the 
distant interactions are considered using delays, but in 
the current study both distant and local interactions are 
assumed to take place simultaneously. 

To realise the long distance rewiring process, a rewiring 
rate r is introduced which denotes the relative rate for 
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FIG. 3: Effective average partial density x'^ (red line 
with crosses) and KMC average partial density {xi) 
(black line with circles) function of the imposed 
reactive reactive coupling pi . The MF partial density 

xi (using the original pi values) is also plotted for 
comparison. Other parameter values are: p2 = 0.016 
and pz = 0.026 and the system size is i ~ 2^. KMC 
averages are taken over 10 runs. 



long distance versus local interactions. For the local 
(short distance) reaction rates the reaction rates are de- 
noted as Pi, i ~ 1,2,3, as before, while for the long dis- 
tance reaction rates the parameters are denoted as Pi^iong , 
i = l,2,3. 

Having chosen the values of r, pi and Pi^iong, « = 1, 2, 3 
the modifications to the KMC algorithm arc straight for- 
ward. With probability 1 — r the classical KMC as de- 
scribed in section IIIBI is realised. With probability r the 
same algorithm holds, but now the neighbours are ran- 
domly chosen between all particles in the system and the 
reactions are realised with rates Pi^iong- In one ETS ei- 
ther a local event or a distant event may take place, thus 
the ratio of number of ETS where local interactions take 
place to the number of ETS where long distance events 
happen is r/(l — r). 

To find out how the rewiring process modulates the 
local process we choose to work with the same reac- 
tion rates, i.e. {pi,P2,Pz) = {Pl. long jP2, long, PS, long) = 

(0.95,0.016,0.026). In this case, while the local interac- 
tions drive the steady state away from the MF solution, 
the rewiring process drives the system towards the MF 
because the "nearest neighbours" are drawn with equal 
probability with the entire system. 

In Fig. m the average position of the xi variable is 
shown as a function of the rewiring rate r. It is first 
noted that the introduction of a rewiring process intro- 
duces a nontrivial deviation on the average values as a 
function of the rate r. The average partial density xi ini- 
tially decreases as the rate of rewiring increases up to a 
critical rate rc ~ 0.8. After this value the average partial 
density xi starts increasing abruptly. The interpretation 
is that for small values of the rewiring rate r the local 
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properties dominate, while for large values of r the long 
distance interactions dominate, stands for the specific 
rewiring value where the global properties take over. The 



ample, for the working parameter set, the minimum of 




FIG. 4: Average densities Xi as a function of the 
rewiring rate r. The black solid line denote the KMC 
simulations results while the red dashed line denotes the 
effective MF values of the species concentration xi. The 
parameter values are: (pi,P2,P3) = (0.95,0.016,0.026) 
and the system size is L = 2^. Averages are taken over 
10 runs and the error bars are designated in the plot. 

In the same plot the average xi densities are computed 
using the effective parameters Pi,eff via Eq. [S] With 
the use of the effective parameters the average popula- 
tion densities are very well predicted by the MF equa- 
tions. This comes as no surprise, after the successful 
results of the previous section, where the effective MF 
theory predicted well the local steady state concentra- 
tions. Here, the introduction of the long distance KMC 
rewiring drives the system towards MF and thus the ef- 
fective MF approximation has a higher advantage in the 
description of the system over the previous, purely local 
KMC approach. 

It it possible to obtain the minimum of curve|3]if we as- 
sume that the overall collective effects due to the stochas- 
tic KMC process and to the long distance reactions is a 
multiplicative factor a, which facilitates the pi rate as 
a'^Pi, the p2 rate as a^p2 and the p^ rate as a^ps- Then 
the critical value Oc could be obtained as the minimum 
with respect to a of Eq. [5] (with the a— modified rates), 
i.e. the minimum with respect to a of 
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derivative of the a;™°'^-component in Eq. [8] with respect 
to a vanishes. In Fig. [5] the derivative is plotted with 
respect to a. It crosses the y = 0-axis at a single value 
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FIG. 5: The derivative dx'{'°''/da as a function of a 



V. KINETIC LATTICE MONTE CARLO WITH 
STOCHASTIC DIFFUSIVE REWIRING 

Long distance, stochastic motion (or long distance dif- 
fusion) may take place in networks with mobile species 
(l9| . As an example, we give the long distance migration 
in ecology, especially in bird motion. In such systems 
species migrate in variable distances and they react lo- 
cally wherever they land. There are many systems in 
which long distance reaction and long distance diffusion 
are both possible. A common one is opinion dynamics 
where the individuals can interact with other distant in- 
dividuals (over the phone or via the Internet) and can 
also relocate and then react locally. Another such net- 
work is the sales network where potential clients can be 
contacted either by Internet (long distance reaction) or 
by salesmen (long distance motion-diffusion). 

In this section we attempt to compare how the long 
distance reaction and long distance diffusion mechanisms 
affect the output of the reactive dynamics. Both mech- 
anisms, long distance reaction and long distance diffu- 
sion, cause a mixing in the system and one could assume 
that the two processes lead to very similar output results. 
Nonetheless, this is not the case and, as it will be shown 
later in this section, the two processes combined with the 
local reactions give very different overall dynamics. 

To compare with the previous section, IIVI we only con- 
sider here local reactions together with long distance dif- 
fusion mechanisms. The local reaction are realised ac- 
cording to the KMC algorithm described in section Hi Bl 
while all the local reaction parameters are the same as 
in sections IIIBI and IIIIl The local reactions are governed 
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by rates pi, while the long distance diffusion rates are 
denoted by Pi.diff- In general the rates Pi^diff may be 
different for the different species. In the current study, 
in order to reduce the number of parameters we assume 
that all species diffuse with the same rate pd- 

Having chosen the values of pd and pi, i = 1,2,3, 
the modifications to the KMC algorithm of section III Bl 
are straight forward. Once a lattice site is selected the 
hosted particle reacts according to the classical KMC as 
described in section Hi Bl with probability 1 — pd- Other- 
wise, with probability pd the hosted particle exchanges 
position with another, randomly selected, particle on the 
lattice. In one ETS either a local reaction event or a 
long distance diffusion event may take place, thus the ra- 
tio of number of ETS where local interactions take place 
to the number of ETS where long distance events happen 
is Pd/(1 — Pd), in direct similarity with the long distance 
reaction process presented in the previous section. Sec. 
IIVI In addition, for comparison with previous results we 
choose to work with the same local reaction rates, i.e. 
{pi,P2,P3) = (0.95,0.016,0.026) and with variable pd- 
The diffusion range Id for the simulations presented here 
are set to Id ~ L/2, but similar results are obtained for 
all Id, which are greater than the size of the local oscil- 
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FIG. 6: Average partial densities xi as a function of the 
long distance reactive coupling pd- The black solid line 

denotes the KMC simulations results, and the red 
dashed curve denotes the Effective MF values of xi. The 
parameter values are {pi,P2,P3) = (0.95,0.016,0.026) 
and the system size is L = 2^. Averages are taken over 
10 configurations and the error bars are shown. 

As can be seen in Fig. [6] there is a decrease in the 
rate of production of Xi particles, which takes place al- 
ready for small values of diffusive mixing, while in the 
case of reactive mixing the decrease is more gradual (sec 
Fig. U). When Pd = ^ only diffusive mixing takes 
place, the effective reaction rates Pi.eff arc and the 
partial densities remain unaltered. That is, for pd = I, 
xi = xi(t = 0), X2 = X2{t = 0), s = s{t = 0). For 



equiprobable initial conditions which arc often used in 
the KMC process, the final states is xi = a;2 = s = 1/3 
for Pd = 1. 

In the same figure the red dashed line represents the 
predictions of the effective MF approach. Unlike in the 
case of reactive mixing the effective MF theory fails to 
predict the resulting rates. This is because the intro- 
duced diffusion is a new process which is not taken into 
account by the effective MF model. And although long 
distance reaction and long distance diffusion both induce 
mixing effects, the details of each processes lead to differ- 
ent outputs in the steady state production of the three 
species. 

VI. CONCLUSIONS 

In this work we study the impact of long distance reac- 
tion and diffusion effects on the steady state of a process 
governed by local Limit Cycle Dynamics. The process 
is realised on a square lattice using Kinetic Monte Carlo 
simulations. 

For the case where only local interactions are consid- 
ered, it is shown that the calculation of a posteriori effec- 
tive reaction rates, allows to get a good approximation 
of the local system steady state properties, using the MF 
effective rates. Further more it is shown, that when long 
distance reactions together with local ones take place, the 
effective MF approximation still faithfully describes the 
steady state properties of the system. 

When long distance diffusion on the local reactive dy- 
namics are added the effective MF description is no more 
valid. This is because the diffusion is a new process which 
is not taken into account by the original MF equations 
which include only reactive terms. Thus the effective 
MF, based only on the reactive terms fail to describe the 
composite system. An extended MF model needs to be 
devised which includes also diffusion effects as a better 
candidate on which to base an effective MF theory to de- 
scribe the dynamics of the LLC model together with the 
long distant diffusion terms. 

In the current study, we have only discussed the prop- 
erties of the steady state and not the dynamics around 
it. In future studies the influence of the long distance 
reactive and diffusive processes on the dynamics of the 
limit cycle, on synchronisation effects and the position of 
the bifurcation point need to be addressed. 
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